Integrating sex-bias into studies of archaic introgression on chromosome X.

Evidence of interbreeding between archaic hominins and humans comes from methods that infer the locations of segments of archaic haplotypes, or 'archaic coverage' using the genomes of people living today. As more estimates of archaic coverage have emerged, it has become clear that most of this coverage is found on the autosomes- very little is retained on chromosome X. Here, we summarize published estimates of archaic coverage on autosomes and chromosome X from extant human samples. We find on average 7 times more archaic coverage on autosomes than chromosome X, and identify broad continental patterns in this ratio: greatest in European samples, and least in South Asian samples. We also perform extensive simulation studies to investigate how the amount of archaic coverage, lengths of coverage, and rates of purging of archaic coverage are affected by sex-bias caused by an unequal sex ratio within the archaic introgressors. Our results generally confirm that, with increasing male sex-bias, less archaic coverage is retained on chromosome X. Ours is the first study to explicitly model such sex-bias and its potential role in creating the dearth of archaic coverage on chromosome X.

High-throughput Oligopaint screen identifies druggable 3D genome regulators.

The human genome functions as a three-dimensional chromatin polymer, driven by a complex collection of chromosome interactions1-3. Although the molecular rules governing these interactions are being quickly elucidated, relatively few proteins regulating this process have been identified. Here, to address this gap, we developed high-throughput DNA or RNA labelling with optimized Oligopaints (HiDRO)-an automated imaging pipeline that enables the quantitative measurement of chromatin interactions in single cells across thousands of samples. By screening the human druggable genome, we identified more than 300 factors that influence genome folding during interphase. Among these, 43 genes were validated as either increasing or decreasing interactions between topologically associating domains. Our findings show that genetic or chemical inhibition of the ubiquitous kinase GSK3A leads to increased long-range chromatin looping interactions in a genome-wide and cohesin-dependent manner. These results demonstrate the importance of GSK3A signalling in nuclear architecture and the use of HiDRO for identifying mechanisms of spatial genome organization.

Mitotic clustering of pulverized chromosomes from micronuclei.

Complex genome rearrangements can be generated by the catastrophic pulverization of missegregated chromosomes trapped within micronuclei through a process known as chromothripsis1-5. As each chromosome contains a single centromere, it remains unclear how acentric fragments derived from shattered chromosomes are inherited between daughter cells during mitosis6. Here we tracked micronucleated chromosomes with live-cell imaging and show that acentric fragments cluster in close spatial proximity throughout mitosis for asymmetric inheritance by a single daughter cell. Mechanistically, the CIP2A-TOPBP1 complex prematurely associates with DNA lesions within ruptured micronuclei during interphase, which poises pulverized chromosomes for clustering upon mitotic entry. Inactivation of CIP2A-TOPBP1 caused acentric fragments to disperse throughout the mitotic cytoplasm, stochastically partition into the nucleus of both daughter cells and aberrantly misaccumulate as cytoplasmic DNA. Mitotic clustering facilitates the reassembly of acentric fragments into rearranged chromosomes lacking the extensive DNA copy-number losses that are characteristic of canonical chromothripsis. Comprehensive analysis of pan-cancer genomes revealed clusters of DNA copy-number-neutral rearrangements-termed balanced chromothripsis-across diverse tumour types resulting in the acquisition of known cancer driver events. Thus, distinct patterns of chromothripsis can be explained by the spatial clustering of pulverized chromosomes from micronuclei.

Recombination between heterologous human acrocentric chromosomes.

The short arms of the human acrocentric chromosomes 13, 14, 15, 21 and 22 (SAACs) share large homologous regions, including ribosomal DNA repeats and extended segmental duplications1,2. Although the resolution of these regions in the first complete assembly of a human genome-the Telomere-to-Telomere Consortium's CHM13 assembly (T2T-CHM13)-provided a model of their homology3, it remained unclear whether these patterns were ancestral or maintained by ongoing recombination exchange. Here we show that acrocentric chromosomes contain pseudo-homologous regions (PHRs) indicative of recombination between non-homologous sequences. Utilizing an all-to-all comparison of the human pangenome from the Human Pangenome Reference Consortium4 (HPRC), we find that contigs from all of the SAACs form a community. A variation graph5 constructed from centromere-spanning acrocentric contigs indicates the presence of regions in which most contigs appear nearly identical between heterologous acrocentric chromosomes in T2T-CHM13. Except on chromosome 15, we observe faster decay of linkage disequilibrium in the pseudo-homologous regions than in the corresponding short and long arms, indicating higher rates of recombination6,7. The pseudo-homologous regions include sequences that have previously been shown to lie at the breakpoint of Robertsonian translocations8, and their arrangement is compatible with crossover in inverted duplications on chromosomes 13, 14 and 21. The ubiquity of signals of recombination between heterologous acrocentric chromosomes seen in the HPRC draft pangenome suggests that these shared sequences form the basis for recurrent Robertsonian translocations, providing sequence and population-based confirmation of hypotheses first developed from cytogenetic studies 50 years ago9.

Whole-genome doubling drives oncogenic loss of chromatin segregation.

Whole-genome doubling (WGD) is a recurrent event in human cancers and it promotes chromosomal instability and acquisition of aneuploidies1-8. However, the three-dimensional organization of chromatin in WGD cells and its contribution to oncogenic phenotypes are currently unknown. Here we show that in p53-deficient cells, WGD induces loss of chromatin segregation (LCS). This event is characterized by reduced segregation between short and long chromosomes, A and B subcompartments and adjacent chromatin domains. LCS is driven by the downregulation of CTCF and H3K9me3 in cells that bypassed activation of the tetraploid checkpoint. Longitudinal analyses revealed that LCS primes genomic regions for subcompartment repositioning in WGD cells. This results in chromatin and epigenetic changes associated with oncogene activation in tumours ensuing from WGD cells. Notably, subcompartment repositioning events were largely independent of chromosomal alterations, which indicates that these were complementary mechanisms contributing to tumour development and progression. Overall, LCS initiates chromatin conformation changes that ultimately result in oncogenic epigenetic and transcriptional modifications, which suggests that chromatin evolution is a hallmark of WGD-driven cancer.

HPV integration generates a cellular super-enhancer which functions as ecDNA to regulate genome-wide transcription.

Human papillomavirus (HPV) integration is a critical step in cervical cancer development; however, the oncogenic mechanism at the genome-wide transcriptional level is still poorly understood. In this study, we employed integrative analysis on multi-omics data of six HPV-positive and three HPV-negative cell lines. Through HPV integration detection, super-enhancer (SE) identification, SE-associated gene expression and extrachromosomal DNA (ecDNA) investigation, we aimed to explore the genome-wide transcriptional influence of HPV integration. We identified seven high-ranking cellular SEs generated by HPV integration in total (the HPV breakpoint-induced cellular SEs, BP-cSEs), leading to intra-chromosomal and inter-chromosomal regulation of chromosomal genes. The pathway analysis revealed that the dysregulated chromosomal genes were correlated to cancer-related pathways. Importantly, we demonstrated that BP-cSEs existed in the HPV-human hybrid ecDNAs, explaining the above transcriptional alterations. Our results suggest that HPV integration generates cellular SEs that function as ecDNA to regulate unconstrained transcription, expanding the tumorigenic mechanism of HPV integration and providing insights for developing new diagnostic and therapeutic strategies.

The p-Arms of Human Acrocentric Chromosomes Play by a Different Set of Rules.

The p-arms of the five human acrocentric chromosomes bear nucleolar organizer regions (NORs) comprising ribosomal gene (rDNA) repeats that are organized in a homogeneous tandem array and transcribed in a telomere-to-centromere direction. Precursor ribosomal RNA transcripts are processed and assembled into ribosomal subunits, the nucleolus being the physical manifestation of this process. I review current understanding of nucleolar chromosome biology and describe current exploration into a role for the NOR chromosomal context. Full DNA sequences for acrocentric p-arms are now emerging, aided by the current revolution in long-read sequencing and genome assembly. Acrocentric p-arms vary from 10.1 to 16.7 Mb, accounting for ∼2.2% of the genome. Bordering rDNA arrays, distal junctions, and proximal junctions are shared among the p-arms, with distal junctions showing evidence of functionality. The remaining p-arm sequences comprise multiple satellite DNA classes and segmental duplications that facilitate recombination between heterologous chromosomes, which is likely also involved in Robertsonian translocations.

Structural variants drive context-dependent oncogene activation in cancer.

Higher-order chromatin structure is important for the regulation of genes by distal regulatory sequences1,2. Structural variants (SVs) that alter three-dimensional (3D) genome organization can lead to enhancer-promoter rewiring and human disease, particularly in the context of cancer3. However, only a small minority of SVs are associated with altered gene expression4,5, and it remains unclear why certain SVs lead to changes in distal gene expression and others do not. To address these questions, we used a combination of genomic profiling and genome engineering to identify sites of recurrent changes in 3D genome structure in cancer and determine the effects of specific rearrangements on oncogene activation. By analysing Hi-C data from 92 cancer cell lines and patient samples, we identified loci affected by recurrent alterations to 3D genome structure, including oncogenes such as MYC, TERT and CCND1. By using CRISPR-Cas9 genome engineering to generate de novo SVs, we show that oncogene activity can be predicted by using 'activity-by-contact' models that consider partner region chromatin contacts and enhancer activity. However, activity-by-contact models are only predictive of specific subsets of genes in the genome, suggesting that different classes of genes engage in distinct modes of regulation by distal regulatory elements. These results indicate that SVs that alter 3D genome organization are widespread in cancer genomes and begin to illustrate predictive rules for the consequences of SVs on oncogene activation.

Dicentric chromosome assay using a deep learning-based automated system.

The dicentric chromosome assay is the "gold standard" in biodosimetry for estimating radiation exposure. However, its large-scale deployment is limited owing to its time-consuming nature and requirement for expert reviewers. Therefore, a recently developed automated system was evaluated for the dicentric chromosome assay. A previously constructed deep learning-based automatic dose-estimation system (DLADES) was used to construct dose curves and calculate estimated doses. Blood samples from two donors were exposed to cobalt-60 gamma rays (0-4 Gy, 0.8 Gy/min). The DLADES efficiently identified monocentric and dicentric chromosomes but showed impaired recognition of complete cells with 46 chromosomes. We estimated the chromosome number of each "Accepted" sample in the DLADES and sorted similar-quality images by removing outliers using the 1.5IQR method. Eleven of the 12 data points followed Poisson distribution. Blind samples were prepared for each dose to verify the accuracy of the estimated dose generated by the curve. The estimated dose was calculated using Merkle's method. The actual dose for each sample was within the 95% confidence limits of the estimated dose. Sorting similar-quality images using chromosome numbers is crucial for the automated dicentric chromosome assay. We successfully constructed a dose-response curve and determined the estimated dose using the DLADES.

Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Imputation-free reconstructions of three-dimensional chromosome architectures in human diploid single-cells using allele-specified contacts.

Single-cell Hi-C analysis of diploid human cells is difficult because of the lack of dense chromosome contact information and the presence of homologous chromosomes with very similar nucleotide sequences. Thus here, we propose a new algorithm to reconstruct the three-dimensional (3D) chromosomal architectures from the Hi-C dataset of single diploid human cells using allele-specific single-nucleotide variations (SNVs). We modified our recurrence plot-based algorithm, which is suitable for the estimation of the 3D chromosome structure from sparse Hi-C datasets, by newly incorporating a function of discriminating SNVs specific to each homologous chromosome. Here, we eventually regard a contact map as a recurrence plot. Importantly, the proposed method does not require any imputation for ambiguous segment information, but could efficiently reconstruct 3D chromosomal structures in single human diploid cells at a 1-Mb resolution. Datasets of segments without allele-specific SNVs, which were considered to be of little value, can also be used to validate the estimated chromosome structure. Introducing an additional mathematical measure called a refinement further improved the resolution to 40-kb or 100-kb. The reconstruction data supported the notion that human chromosomes form chromosomal territories and take fractal structures where the dimension for the underlying chromosome structure is a non-integer value.

Short arms of human acrocentric chromosomes and the completion of the human genome sequence.

The complete, ungapped sequence of the short arms of human acrocentric chromosomes (SAACs) is still unknown almost 20 years after the near completion of the Human Genome Project. Yet these short arms of Chromosomes 13, 14, 15, 21, and 22 contain the ribosomal DNA (rDNA) genes, which are of paramount importance for human biology. The sequences of SAACs show an extensive variation in the copy number of the various repetitive elements, the full extent of which is currently unknown. In addition, the full spectrum of repeated sequences, their organization, and the low copy number functional elements are also unknown. The Telomere-to-Telomere (T2T) Project using mainly long-read sequence technology has recently completed the assembly of the genome from a hydatidiform mole, CHM13, and has thus established a baseline reference for further studies on the organization, variation, functional annotation, and impact in human disorders of all the previously unknown genomic segments, including the SAACs. The publication of the initial results of the T2T Project will update and improve the reference genome for a better understanding of the evolution and function of the human genome.

The complete sequence of a human genome.

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion-base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

Most complete human genome yet is revealed.

Unusual cell line helps sequencers read previously indecipherable stretches of DNA.

The genomic and spatial mobility of extrachromosomal DNA and its implications for cancer therapy.

Extrachromosomal DNA (ecDNA) amplification has been observed in at least 30 different cancer types and is associated with worse patient outcomes. This has been linked to increased oncogene dosage because both oncogenes and associated enhancers can occupy ecDNA. New data challenge the view that only oncogene dosage is affected by ecDNA, and raises the possibility that ecDNA could disrupt genome-wide gene expression. Recent investigations suggest that ecDNA localizes to specialized nuclear bodies (hubs) in which they can act in trans as ectopic enhancers for genes on other ecDNA or chromosomes. Moreover, ecDNA can reintegrate into the genome, possibly further disrupting the gene regulatory landscape in tumor cells. In this Perspective, we discuss the emerging properties of ecDNA and highlight promising avenues to exploit this new knowledge for the development of ecDNA-directed therapies for cancer.

Non-Random Pattern of Integration for Epstein-Barr Virus with Preference for Gene-Poor Genomic Chromosomal Regions into the Genome of Burkitt Lymphoma Cell Lines.

BACKGROUND: Epstein-Barr virus (EBV) is an oncogenic virus found in about 95% of endemic Burkitt lymphoma (BL) cases. In latently infected cells, EBV DNA is mostly maintained in episomal form, but it can also be integrated into the host genome, or both forms can coexist in the infected cells. METHODS: In this study, we mapped the chromosomal integration sites of EBV (EBV-IS) into the genome of 21 EBV+ BL cell lines (BL-CL) using metaphase fluorescence in situ hybridization (FISH). The data were used to investigate the EBV-IS distribution pattern in BL-CL, its relation to the genome instability, and to assess its association to common fragile sites and episomes. RESULTS: We detected a total of 459 EBV-IS integrated into multiple genome localizations with a preference for gene-poor chromosomes. We did not observe any preferential affinity of EBV to integrate into common and rare fragile sites or enrichment of EBV-IS at the chromosomal breakpoints of the BL-CL analyzed here, as other DNA viruses do. CONCLUSIONS: We identified a non-random integration pattern into 13 cytobands, of which eight overlap with the EBV-IS in EBV-transformed lymphoblastoid cell lines and with a preference for gene- and CpGs-poor G-positive cytobands. Moreover, it has been demonstrated that the episomal form of EBV interacts in a non-random manner with gene-poor and AT-rich regions in EBV+ cell lines, which may explain the observed affinity for G-positive cytobands in the EBV integration process. Our results provide new insights into the patterns of EBV integration in BL-CL at the chromosomal level, revealing an unexpected connection between the episomal and integrated forms of EBV.

Mutations along human chromosomes: How randomly scattered are they?

The diversity of mutations in human chromosomes is nowadays very well documented. The mutations characterize populations in the world as well as genetic causes of diseases. In the approach that we follow, we study the patterns of gaps between mutations by means of the rescaled range analysis and the fractal dimension estimates. The results for chromosomes 1 to 22 and X indicate the existence of the so-called Hurst phenomenon in all of them. The interpretation of this outcome entails the presence of long-range correlations and we propose an explanation based on the genomic feature dubbed linkage disequilibrium, a nonrandom association of alleles at different loci. An unexpected outcome is the noteworthy uniform reduction in the Hurst phenomenon when considering the centimorgan metric instead of base position units. By contrast, such uniform reduction is not observed with the fractal dimension values.

cGAS guards against chromosome end-to-end fusions during mitosis and facilitates replicative senescence.

As a sensor of cytosolic DNA, the role of cyclic GMP-AMP synthase (cGAS) in innate immune response is well established, yet how its functions in different biological conditions remain to be elucidated. Here, we identify cGAS as an essential regulator in inhibiting mitotic DNA double-strand break (DSB) repair and protecting short telomeres from end-to-end fusion independent of the canonical cGAS-STING pathway. cGAS associates with telomeric/subtelomeric DNA during mitosis when TRF1/TRF2/POT1 are deficient on telomeres. Depletion of cGAS leads to mitotic chromosome end-to-end fusions predominantly occurring between short telomeres. Mechanistically, cGAS interacts with CDK1 and positions them to chromosome ends. Thus, CDK1 inhibits mitotic non-homologous end joining (NHEJ) by blocking the recruitment of RNF8. cGAS-deficient human primary cells are defective in entering replicative senescence and display chromosome end-to-end fusions, genome instability and prolonged growth arrest. Altogether, cGAS safeguards genome stability by controlling mitotic DSB repair to inhibit mitotic chromosome end-to-end fusions, thus facilitating replicative senescence.

Clinical utility of chromosomal microarray in establishing clonality and high risk features in patients with Richter transformation.

Richter transformation (RT) refers to the development of an aggressive lymphoma in patients with pre-existing chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL/SLL). It carries a poor prognosis secondary to poor response to therapy or rapid disease relapse. Currently there are no randomized trials to guide treatment. Therapeutic decisions are often influenced by the presence or absence of a clonal relationship between the underlying CLL/SLL and the new lymphoma given the poor prognosis of patients with clonally related RT. Chromosomal microarray analysis (CMA) can help to establish clonality while also detecting genomic complexity and clinically relevant genetic variants such as loss of CDKN2A and/or TP53. As a result, CMA has potential prognostic and therapeutic implications. For this study, CMA results from patients with Richter transformation were evaluated in paired CLL/SLL and transformed lymphoma samples. CMA revealed that 86% of patients had common aberrations in the two samples indicating evidence of common clonality. CMA was also useful in detecting aberrations associated with a poor prognosis in 71% of patients with RT. This study highlights the potential clinical utility of CMA to investigate the clonal relationship between CLL/SLL and RT, provide prognostic information, and possibly guide therapeutic decision making for patients with Richter transformation.

CIR-Net: Automatic Classification of Human Chromosome Based on Inception-ResNet Architecture.

BACKGROUND: In medicine, karyotyping chromosomes is important for medical diagnostics, drug development, and biomedical research. Unfortunately, chromosome karyotyping is usually done by skilled cytologists manually, which requires experience, domain expertise, and considerable manual efforts. Therefore, automating the karyotyping process is a significant and meaningful task. METHOD: This paper focuses on chromosome classification because it is critical for chromosome karyotyping. In recent years, deep learning-based methods are the most promising methods for solving the tasks of chromosome classification. Although the deep learning-based Inception architecture has yielded state-of-the-art performance in the 2015 ILSVRC challenge, it has not been used in chromosome classification tasks so far. Therefore, we develop an automatic chromosome classification approach named CIR-Net based on Inception-ResNet which is an optimized version of Inception. However, the classification performance of origin Inception-ResNet on the insufficient chromosome dataset still has a lot of capacity for improvement. Further, we propose a simple but effective augmentation method called CDA for improving the performance of CIR-Net. RESULTS: The experimental results show that our proposed method achieves 95.98 percent classification accuracy on the clinical G-band chromosome dataset whose training dataset is insufficient. Moreover, the proposed augmentation method CDA improves more than 8.5 percent (from 87.46 to 95.98 percent) classification accuracy comparing to other methods. In this paper, the experimental results demonstrate that our proposed method is recent the most effective solution for solving clinical chromosome classification problems in chromosome auto-karyotyping on the condition of the insufficient training dataset. Code and Dataset are available at https://github.com/CloudDataLab/CIR-Net.